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ABSTRACT 

We have performed a detailed systematic search for multiperiodicity in the Population I Cepheids 
of the Large Magellanic Cloud. In this process we have identified for the first time several new types 
of Cepheid pulsational behaviour. We have found two triple-mode Cepheids pulsating simultaneously 
in the first three radial overtones. In 9% of the first overtone Cepheids we have detected weak, but well 
resolved secondary periodicities. They appear either very close to the primary pulsation frequency 
or at a much higher frequency with a characteristic period ratio of 0.60 — 0.64 . In either case, the 
secondary periodicities must correspond to nonradial modes of oscillation. This result presents a 
major challenge to the theory of stellar pulsations, which predicts that such modes should not be 
exited in Cepheid variables. Nonradial modes have also been found in three of the fundamental/first 
overtone double-mode Cepheids, but no such oscillations have been detected in single mode Cepheids 
pulsating in the fundamental mode. 

In 19% of double-mode Cepheids pulsating in the first two radial overtones (FO/SO type) we 
have detected a Blazhko-type periodic modulation of amplitudes and phases. Both modes are modu- 
lated with a common period, which is always longer than 700 days. Variations of the two amplitudes 
are anticorrelated and maximum of one amplitude always coincides with minimum of the other We 
have compared observations of modulated FO/SO Cepheids with predictions of theoretical models 
of the Blazhko effect, showing that currently most popular models cannot account for properties of 
these stars. We propose that Blazhko effect in FO/SO Cepheids can be explained by a nonstationary 
resonant interaction of one of the radial modes with another, perhaps nonradial, mode of oscillations. 
Key words: Stars: variables: Ceplieids — stars: oscillations: nonradial - stars: oscillations: multi- 
mode - stars: oscillations: Blazhko ejfect — galaxies: Magellanic Clouds 



1. Introduction 



Classical Cepheids and RR Lyrae stars are among the best studied pulsating 
variables. For decades they have been considered primary example of purely ra- 
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dial pulsators, displaying large amplitude oscillations, with only one or two radial 
modes excited. This simple picture was challenged in the end of the last cen- 
tury when Olech et al. (1999) discovered low amplitude nonradial modes in three 
RR Lyr stars of globular cluster M55. In following years it was shown that presence 
of secondary periodicities identified with nonradial modes of oscillations is a very 
common property of RR Lyr variables. They were detected in a substantial fraction 
of stars in every studied RR Lyr population: in LMC (Alcock et al. 2000, 2003; 
Soszyriski et al. 2003; Nagy & Kovacs 2006), SMC (Soszyriski et al. 2002), Galac- 
tic Bulge (Moskalik & Poretti 2002, 2003; Mizerski 2003; CoUinge et al. 2006), 
Galactic field (Wils et al. 2006; Szczygiel & Fabrycky 2007, see also Smith 1995) 
and most recently in CO Cen (Olech & Moskahk 2008). Depending on the stellar 
system, nonradial modes are excited in between 5% and 30% of all RR Lyr stars. 

Motivated by these findings we decided to search for nonradial modes in classi- 
cal (Population I) Cepheids. As far as pulsation properties are concerned, Cepheids 
are very close siblings of RR Lyr stars. It is therefore very interesting to check 
if nonradial modes are excited in these pulsators as well. With this goal in mind, 
we conducted a systematic search for multiperiodicity among Pop. I Cepheids of 
the Large Magellanic Cloud (LMC). The primary source of data for this analy- 
sis was the OGLE-II photometry (Zebrun et al. 2001). The data, collected du- 
ring microlensing survey, covers the timebase of 1000-1200 days, with 250-500 
I-band measurements per star. Such a long, uniform and high quality dataset is 
particularly well suited for our purpose. Partial results of this survey were pub- 
lished by Moskalik, Kolaczkowski & Mizerski (2004, 2006) and by Moskalik & 
Kolaczkowki (2008a,b). In this paper we present a full discussion of our findings 
and provide a complete inventory of newly identified multiperiodic Cepheids in the 
LMC. 

2. Time series analysis of Cepheid lightcurves 

We analyzed all fundamental mode (FU), first overtone (FO) and double-mode 
(FU/FO and FO/SO) LMC Cepheids Usted in the OGLE-II catalogs (Udalski et al. 
1999b; Soszyriski et al. 2000), nearly 1300 stars in total. We omitted all variables 
marked by the authors as FA or BR. The former are most likely Pop. n Cepheids, 
while the status of the latter group is somewhat unclear. We used the I-band OGLE- 
II data obtained with the Difference Image Analysis (or DIA) method. The search 
for secondary periodicities was conducted with a standard consecutive prewhiten- 
ing technique. First, we fitted the data with a Fourier sum describing pulsations 
with the dominant mode: 

lit) = (/) +£A^sin(2ji/:/if + (^^). (1) 
k 

The frequency of the mode, f\ , was also optimized. For double-mode Cepheids, 
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we fitted the data with double-frequency Fourier sum, representing pulsations in 
two dominant (radial) modes: 

I{t) = {!) + T^^jk sin[27i(y/i + kf2)t + ^jk] . (2) 

The residuals of the fit were then searched for secondary frequencies. This was 
done with the Fourier transform, calculated over the range of — 6 c/d. If any new 
periodicity was detected, a new Fourier sum, including all frequencies identified so 
far, was fitted to the data and the fit residuals were examined again. The process 
was stopped when no new significant frequencies appeared. At this stage we per- 
formed data clipping, rejecting all measurements deviating from the fitted function 
by more than 5a, where a is the standard deviation of the fit residuals. OGLE-II 
photometry is very clean and only for a small number of stars 1-3 datapoints were 
removed by this criterion. After data clipping, the frequency analysis was repeated. 

As a final step of analysis we checked all Cepheids with detected secondary 
peaks for a possible lightcurve contamination. It is a well known shortcoming of 
the DIA reduction method, that the measured flux of a star can be contaminated 
by light coming from neighbouring stars {i.e. Mizerski & Bejger 2002; Hartman 
et at. 2004). As a result, some of the discovered secondary frequencies might origi- 
nate not in the Cepheid, but in another variable star in Cepheid's neighbourhood. In 
order to remove such spurious detections we proceeded in the following way. First, 
we identified all variable stars from OGLE-II catalog (Zebruri et al. 2001), which 
are located within r = 1 arcmin of the studied Cepheid. For each of these stars 
we performed frequency analysis and then checked if detected frequencies or their 
aliases appear as secondary peaks in Cepheid's power spectrum. With this proce- 
dure we weeded out 14 contaminated Cepheids. In LMC fields SCI to SC12 we ex- 
tended the search radius to r = 1 .5 arcmin, but no new cases of contamination were 
found. The OGLE-II catalog of variable stars is not complete, though. It misses 
most small amplitude variables, which can contaminate Ughtcurves of nearby stars 
as well. Therefore, we performed additional contamination check, calculating fre- 
quency spectra of all stars (not only known variables) within r = 20 arcsec of the 
Cepheid. No more contaminated Cepheids were rejected at this step. 

3. New double-mode Cepheids 

In the course of our analysis we identified five new "canonical" double-mode 
Cepheids, pulsating in two radial modes. In four of them the fundamental mode 
and the first overtone are excited (FU/FO type), while in the fifth Cepheid the first 
and second overtones are excited (FO/SO type). Basic properties of these variables 
are listed in Table 1. Periods and ampUtudes given in the table (as well as in all 
other tables in this paper) are determined through the least square fits of appropriate 
multifrequency solutions to the data. 
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Table 1 

New double-mode Cepheids in OGLE-II LMC sample 



OGLE ID 


Po [day] 


Pi [day] 


Pi/Po 


AiMo 


Type 


SC2-263415 
SC 11-250925 
SC17-126402'' 
SC20~100652 


3.420377 
7.864120 
3.209709 
3.747624 


2.455273 
5.565078 
2.294942 
2.715290 


0.71784 
0.70765 
0.71500 
0.72454 


21. n 
11.11 

24.31 
16.65 


FU/FO 
FU/FO 
FU/FO 
FU/FO 




Pi [day] 


Pi [day] 


Pl/Pl 






SC20-1 12788* 


0.737740 


0.594287 


0.80555 


0.09 


FO/SO 



" Additional mode at P=l. 429530 day detected (see Section 6). 
* Both radial modes modulated (see Section 7). 



In Fig. 1 and Fig. 2 we present the period ratio vs. period diagrams (so called 
Petersen diagrams) for all double-mode Cepheids detected in the LMC (Alcock 
et al. 1995, 1999, 2003; Soszyriski et al. 2000). For comparison, we also display 
the FU/FO Cepheids observed in the Small IVIagellanic Cloud (Alcock et al. 1997; 
Udalski et al. 1999a) and in the Galaxy (Pardo & Poretti 1997; Antipin 1997a, 
1998; Berdnikov & Turner 1998; Wils & Otero 200413. Because of different metal- 
licities, the P\_/Pq period ratios in these three stellar systems are somewhat dif- 
ferent. Four of the newly identified double-mode Cepheids fit to their respective 
Petersen diagrams very well. For the fifth object, SC20-100652, the period ratio 
Pi/Po is rather high and is close to the value expected for an SMC star. This sug- 
gests that SC20-100652 has a lower metallicity than other LMC Cepheids. We 
also note that, the new FU/FO double-mode pulsator SC 11-250925 has the longest 
periods among currently known variables of this class. 

4. Triple-mode Cepheids 

In two of the previously known FO/SO double-mode Cepheids a third strong 
periodicity was found. The period ratio of P-i/P2 = 0.840, the same in both Ce- 
pheids, identifies the new mode as a third radial overtone. In Fig. 3 we present the 
prewhitening sequence for one of these stars. The third overtone is detected very 
securely, in fact it is stronger than the second overtone. 

Pulsation properties of the two triple-mode Cepheids are listed in Table 2. In 
Fig. 4 we display both variables in the Cepheid P — L plot constructed for extinction 
insensitive Wessenheit index, W/ = /— 1.55 x (V — /) (Madore & Freedman 1991). 
The stars fit the P — L relation very well, proving that apart from multiperiodic 



For DZ CMa pulsation periods have been redetermined with new data of Berdnikov (2008). 
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Table2 

Triple-mode Cepheids in OGLE-II LMC sample 



OGLE ID 


P\ [day] 


Pi [day] 


Pi [day] 


Pi/Px 


P^/Pi 


A2/A1 


A3/A1 


SC3-360128 
SC5-338399 


0.541279 
0.579510 


0.436046 
0.466624 


0.366300 
0.392116 


0.80559 
0.80520 


0.84005 
0.84033 


0.167 
0.111 


0.311 
0.178 



pulsations, they are both normal Pop. I Cepheids belonging to the LMC. 

Discovery of the triple-mode pulsators listed in Table 2 was originally an- 
nounced by Moskalik, Kolaczkowski & Mizerski (2004). Three more triple-mode 
Cepheids have been identified very recently by Soszyhski et al. (2008a), bringing 
the number of objects in the class to five. This new type of multimode Cepheid 
pulsations, although extremely rare, is very important for testing stellar models. 
With three radial modes observed and their periods accurately measured, seismo- 
logical analysis yields tight constraints on Cepheid parameters (mass, luminosity, 
metalhcity). This imposes constraints on Cepheid evolutionary tracks, including 
interesting limits on convective overshooting from the core (MoskaUk & Dziem- 
bowski 2005). 

5. Analysis of FO and FU Cepheids 

5.1. Nonradial modes in first overtone Cepheids 

The OGLE-II catalog of LMC Cepheids (Udalski et al. 1999b) lists 462 first 
overtone variables. We detected resolved secondary frequencies in 42 of them. 
This constitutes 9% of the sample. We consider two frequencies to be resolved, 
if they differ by more than A/ = 2/r, where T is the length of the data. In case 
of OGLE-II photometry, this means A/ > 0.0017 - 0.0020 c/d, depending on the 
star. Our criterion is more conservative than usually adopted in other studies (i.e. 
Alcock et al. 2003; Nagy & Kovacs 2006). 

The complete inventory of resolved FO Cepheids is presented in Table 3. Fol- 
lowing notation originally introduced for the RR Lyr variables (Alcock et al. 2000), 
we call these stars FO-V Cepheids. Consecutive columns of Table 3 give OGLE 
number of the star, primary and secondary periods Pi and Py , frequency difference 
A/ = /v — /i , period ratio P^/Pi and amplitude ratio Ay/Aj . 

In most of the FO-V Cepheids only one secondary frequency was detected, but 
in several variables two frequencies were found. In all cases they have extremely 
low amphtudes. With the exception of a singe star (SC18-208875), the amplitude 
ratio Ay/Ai is always below O.I, with the average value of 0.048. We note, that 
secondary peaks detected in LMC first overtone RR Lyr stars are several times 
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Table3 

FO-v Cepheids in OGLE-H LMC sample 



OGLE ID 



Pi [day] Pv[day] A/[c/d] Pm/Pi Ay/Ai Remarks 



SCl-201683 2.319211 



SCl-229541 

SC2-158672 
SC2-208897 
SC2-283723 
SC3-274410 

SC3-421512 
SC4-36200 

SC4-131738 

SC4-295932 

SC5-75989 

SC5-138031 

SC6-135695 

SC6-135716 

SC6-267289 

SC6-363194 

SC7-344559 

SC8-205108 

SC8-224964 

SC9-2 16934 

SC9-230584 

SClO-95827 
SC10~132645 
SClO-259946 
SC13-165223 



1.216360 

2.313334 
2.417847 
1.308907 
2.446161 

3.186057 
3.476962 

3.122245 
4.166232 
2.238255 
2.681369 
1.864922 
2.838799 
1.845746 
2.797085 
2.062992 
3.515988 
2.866473 
4.001017 
4.395271 

1.974513 
1.576676 
5.075153 
2.043172 



SCI 3-242700 3.452472 



SC14-46315 

SC15-170744 

SC16-37119 

SC16-177823 

SC 16-194279 

SCI 6-230207 

SC17-39484 



2.315774 
4.992984 
2.795951 
1.918044 
2.065767 
3.418163 
3.541888 



SC17-39517 2.113535 



2.340033 
3.069683 
1.359217 
1.477135 
2.346886 
1.425531 
2.981644 
3.099460 
4.652311 
3.912773 
4.057330 
3.780949 
2.633729 
2.949657 
3.249411 
1.949614 
2.780248 
1.906899 
3.727012 
2.140660 
3.905659 
3.189296 
4.997407 
5.672324 
5.860579 
2.070361 
1.594462 
6.183343 
2.039852 
1.230847 
3.262793 
3.363016 
1.442102 
3.148984 
3.520911 
1.892260 
2.146254 
4.416698 
4.063954 
4.104648 
2.149961 
1.326966 



-0.00384 
-0.10542 
-0.08641 
0.24471 
0.01251 
-0.06250 
-0.07342 
-0.08617 
-0.09892 
-0.03203 
-0.04114 
-0.05580 
0.13967 
-0.10775 
-0.06520 
-0.02329 
0.00742 
-0.01738 
-0.08920 
-0.01759 
-0.02884 
-0.03531 
-0.04983 
-0.05122 
-0.05689 
-0.02345 
-0.00708 
-0.03531 
0.00080 
0.32301 
0.01684 
0.00771 
0.26161 
0.11728 
-0.07364 
0.00710 
-0.01815 
-0.06614 
-0.03627 
-0.03871 
-0.00802 
0.28046 



1.00898 
1.32359 
1.11745 
0.63853 
0.97065 
1.08910 
1.21891 
1.26707 
1.46020 
1.12534 
1.16692 
1.21097 
0.63216 
1.31784 
1.21185 
1.04541 
0.97938 
1.03313 
1.33246 
1.03765 
1.11083 
1.11262 
1.24903 
1.29055 
1.33338 
1.04854 
1.01128 
1.21836 
0.99838 
0.60242 
0.94506 
0.97409 
0.62273 
0.63068 
1.25929 
0.98656 
1.03896 
1.29213 
1.14740 
1.15889 
1.01725 
0.62784 



0.048 
0.022 
0.043 
0.032 
0.074 
0.070 
0.043 
0.033 
0.055 
0.054 
0.038 
0.032 
0.023 
0.091 
0.029 
0.037 
0.064 
0.039 
0.043 
0.068 
0.026 
0.041 
0.031 
0.059 
0.028 
0.076 
0.053 
0.039 
0.040 
0.041 
0.060 
0.077 
0.023 
0.055 
0.051 
0.046 
0.059 
0.039 
0.055 
0.045 
0.065 
0.046 
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Table3 
Concluded 



OGLE ID Pi [day] A, [day] A/[c/d] Py/Pi AyjAy Remarks 



SC17-80220 


2.258854 


2.829748 


-0.08931 


1.25274 


0.058 




SC17-146711 


1.883131 


1.137820 


0.34784 


0.60422 


0.049 




SC17-171481 


2.537651 


3.707561 


-0.12435 


1.46102 


0.066 




SC17-211310 


2.023579 


2.087774 


-0.01519 


1.03172 


0.054 


s3 






2.593589 


-0.10861 


1.28168 


0.037 


a 


SC18-144653 


4.158193 


5.222752 


-0.04902 


1.25602 


0.049 




SCI 8-208875 


3.928946 


3.508220 


0.03052 


0.89292 


0.202 




SC19-74265 


2.364944 


2.790232 


-0.06445 


1.17983 


0.046 


a 






3.080875 


-0.09826 


1.30273 


0.060 




SC20-83423 


2.068067 


2.108720 


-0.00932 


1.01966 


0.053 





Remarks: a) Marginal detection; b) Unresolved residual power at primary fre- 
quency; c) Unresolved residual power at secondary frequency; d) Analyzed with 
MACHO data; si) Same as SC6^22324; s2) Same as SC17-33268; s3) Same 
as SC18-29237. 



higlier, witli Ay/A\ = 0.31 on average (Nagy & Kovacs 2006). 

The secondary frequencies in FO-v Cepheids come in two different flavours. 
In 37 variables tliey are located close to the primary (radial) frequency f\ , within 
|A/| < 0.13 c/d. Representative examples of such behaviour are displayed in Fig. 5. 
In 84% of cases, secondary frequencies are lower than the primary one. When two 
secondary peaks are detected, they always appear on the same side of the primary 
pealc. The observed frequency pattern cannot be explained by amplitude or phase 
modulation of the radial mode. A periodic modulation always produces an equally 
spaced frequency multiplet (triplet, quintuplet, etc.), centered on the primary peak. 
It is also easy to check, that measured period ratios are incompatible with those of 
the radial modes. This implies, that secondary frequencies detected close to the ra- 
dial mode of the FO Cepheids must correspond to nonradial modes of oscillations. 
This is the first detection of nonradial modes in this type of stars. 

In seven FO-v Cepheids secondary periodicities of a different type are found: 
we see a high frequency mode, with characteristic period ratio of Py/Pi = 0.60 — 
0.64. Examples of this behaviour are displayed in Fig. 6. The two types of se- 
condary modes are not mutually exclusive. In two objects (SC13-165223 and 
SC 17-395 17) both a high frequency secondary peak and a secondary peak close 
to the primary frequency are present. 

The mysterious period ratio of 0.60 — 0.64 cannot be explained in terms of ra- 
dial modes, neither. Observationally determined period ratio of the third and first 
radial overtones is P3/P1 = 0.6766 - 0.6773 (see Table 2; Soszynski et al. 2008a). 
The value of P4/P1 can be estimated theoretically. Assuming validity of p-mode 
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asymptotic relation A -/a = h-fi,^^ get P^/P^ = {l-P^/Pi)-^ . With P3/P2 = 
0.840 (see Table 2) we find A/^s = 0.862, thus P^/Pi = 0.5835. This estimate 
agrees nicely with Pa/P\ = 0.5898 derived from tentative detection of the fourth 
overtone in one of the LMC Cepheids (Soszynski et al. 2008a). With Py/^i = 
0.60 — 0.64, the high frequency secondary peaks in FO-v Cepheids cannot corre- 
spond to neither the third nor the fourth radial overtone. They have to be attributed 
to nonradial modes. These modes appear at frequencies 0.02 — 0.06 c/d below fre- 
quency of the (unobserved) fourth radial overtone. Interestingly, this frequency 
difference is comparable to A/ for the nonradial modes detected in the vicinity of 
the first overtone. 

In Fig. 7 we show the distribution of frequency differences. A/, for the LMC 
FO-v Cepheids and, for comparison, for the LMC first overtone Blazhko RR Lyr 
stars (RRc-BL stars). The two distributions are very similar. FO-V Cepheids 
show stronger preference for negative A/, but otherwise, in both types of over- 
tone pulsators nonradial modes are found in similar distances from the radial mode. 
The only difference between the two histograms is presence of high frequency se- 
condary modes {Pv/Pi = 0.60 — 0.64), which were detected in FO Cepheids, but 
not in RRc stars. 

5.2. First overtone Cepheids with unresolved residual power 

In 23 first overtone Cepheids we found after prewhitening a significant resi- 
dual power unresolved from primary pulsation frequency. This is a signature of 
a slow amplitude and/or phase variability, occuring on a timescale comparable or 
longer than the ~ 1100 days length of OGLE-II data. We reanalyzed all these stars 
with the MACHO b photometry (AUsman & Axelrod 2001), which offered a much 
longer timebase of 2700-2800 days. MACHO data was cleaned before analysis. We 
removed all datapoints with formal errors more than 5 times larger than the average 
formal error in the dataset. We also rejected all measurements differing from the 
mean brightness by more than 5 standard deviations (of the unphased data). After 
cleaning, the data was analyzed in a standard way. Remnant power was resolved 
with MACHO photometry only in one additional Cepheid, SCI 3-1 65223. This 
object is included in Table 3. Remaining 22 Cepheids with unresolved residual 
power are listed in Table 4. Following notation introduced by Alcock et al. (2000), 
we call these variables FO-PC Cepheids. 

For closer examination of slow amplitude/phase variability in FO-PC Cepheids 
we divided the data into 10 subsets, each spanning 10% of the total timebase, and 
then fitted Eq.(l) to each subset separately. This way we were able to follow the 
time evolution of amplitude and phase of the radial mode. In 8 stars only parabolic 
trends in pulsation phases were found. These Cepheids most likely undergo a se- 
cular period change. In remaining 14 FO-PC Cepheids variabihty of amplitudes 
is clearly seen, in most cases associated with variability of phases. In Fig. 8 we 
show typical examples of both behaviours. With currently available data we caimot 
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Table4 

FO-PC Cepheids in OGLE-H LMC sample 



OGLE ID Pi [day] Ai [mag] Remarks 



SCl-158032 


2.161982 


— 


ab,c 


SCl-324986 


2.209283 


0.173 


ab,sl 


SC2-248540 


1.823789 


— 


a,d 


SC4~152290 


1.887862 


0.082 


ab 


SC8-76179 


2.462125 


0.183 


b,s2 


SC9-153398 


0.947292 


0.087 


ab 


SC 10-269402 


1.912859 


0.114 


ab 


SCI 1-338308 


1.743443 


0.053 


ab 


SCI 3-74204 


1.369660 


0.079 


a 


SC14-148402 


2.193230 


0.209 


b 


SC15-31558 


1.554964 


0.046 


a 


SC 15-158999 


1.100853 


0.266 


b 


SC16-99240 


2.437482 


0.152 


b 


SC17-15475 


2.012359 


0.200 


ab 


SC 17-1 17748 


3.572147 


0.160 


b 


SC17-122399 


2.783800 


0.067 


a 


SC17-214859 


1.954504 


0.206 


b,s3 


SC18-141019 


2.270553 


0.192 


b 


SCI 8-202349 


2.432908 


0.179 


b,s4 


SC2 1-85282 


3.372482 


0.167 


ab 


SC21-106503 


2.535978 


0.184 


ab 


SC21-136158 


3.774050 


0.153 


ab 



Remarks: a) Amplitude variability; b) Phase 
variability; c) no residual power detected in 
MACHO data; d) no MACHO data available; 
si) Same as SC16-57446; s2) Same as SC9- 
372261; s3) Same as SC18-33591; s4) Same 
as SC 19^8662. 

malce any firm statement about tbie cause of ttie observed long term trends. We only 
note, that in all FO-PC Cepheids the observed changes of phases and ampUtudes 
are orders of magnitude too fast to be explained by stellar evolution. We recall, 
that fast phase variations of nonevolutionary origin are not unique to Cepheids of 
the LMC, they are also observed in many FO Cepheids in the Galaxy (Berdnikov 
et al. 1997, their Fig. 3). On the other hand, secular amplitude variations are rare in 
Galactic FO Cepheids, with Polaris being the only well documented case {e.g. Lee 
et al. 2008; Bruntt et al. 2008). 

5.3. Fundamental mode Cepheids 

The OGLE-n catalog of LIVIC Cepheids (Udalski et al. 1999b) list 718 funda- 
mental mode pulsators. We searched all of them for secondary periodicities. We 
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found no nonradial modes in the FU Cepheids of the LMC. This result is highly 
significant, considering that the FU Cepheid sample is much larger than the FO 
Cepheid sample. Long-term amplitude or phase variability is also extremely rare 
in the FU Cepheids. We found this behaviour only in one star (SC17-94969). 

5.4. Incidence rate of nonradial modes 

In Table 5 we present the inventory of different types of pulsators identified in 
our survey of FO and FU Cepheids in the LMC. For each type of variability we give 
the number of stars found, A'^, , and the estimated incidence rate in the population 
with its statistical error. The errors are calculated from the assumption of Poisson 
distribution {cf. Alcock et al. 2003): 

^^.■/^ = ^V — — ' 

where A'^ is the total number of FO or FU Cepheids. 

Tables 

Variability types in FO and FU Cepheids of OGLE-II LMC sample 



Type 


Description 


Number 


Inc. rate 


FO-S 


Single-periodic overtone Cepheid 


397 


85.9 ±1.6% 


FO-v 


FO with nonradial modes 


42 


9.1 ±1.3% 


FO-PC 


FO with variable phase/amplitude 


22 


4.8 ±1.0% 


FO-Misc 


FO with some miscellany 


1 


0.2 ±0.2% 


FU-S 


Single-periodic fundamental Cepheid 


716 


99.6 ±0.2% 


FU-V 


FU with nonradial modes 






FU-PC 


FU with variable phase/amplitude 


1 


0.1 ±0.1% 


FXJ-Misc 


FU with some miscellany 


1 


0.1 ±0.1% 



The most significant result of our survey is the detection of nonradial modes 
in LMC Cepheids, but only in those which pulsate in the first overtone. In Fig. 9 
we display all Cepheids with nonradial modes on the Wessenheit index P — L dia- 
gram. All variables classified as FO-V Cepheids are indeed firmly located on the 
first overtone sequence. Their disttibution with the pulsation period is not uniform, 
though. FO-V pulsators are not found for Pi < 1.2 day. This suggests that the inci- 
dence rate of nonradial modes might depend on the period. To test this hypothesis 
we divided the FO Cepheid population into several period bins and estimated the in- 
cidence rates in each bin separately. Results are displayed in Fig. 10. The incidence 
rate of FO-v pulsators systematically increases with period, reaching 19 ±5% for 
log Pi > 0.5. We think, that this is most likely a selection effect. The nonradial 
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modes, which have very small amplitudes, are progressively more difficult to de- 
tect in Cepheids with shorter periods i.e. with lower luminosities. This interpreta- 
tion imphes, that the true incidence rate of FO-v Cepheids might be significantly 
higher than 9% given in Table 5. 

5.5. Comparison with RR Lyrae stars 

In the last decade nonradial modes were detected in many RR Lyr stars belong- 
ing to various stellar systems. Since pulsations of Cepheids and of RR Lyr stars 
are in many ways similar, it is interesting to compare the properties of nonradial 
modes in these two types of variables. One similarity has already been discussed 
in Section 5.1: in both types of pulsators, nonradial modes appear in essentially the 
same frequency distances from the radial mode (Fig. 7). However, there are also 
striking differences between the two groups of stars: 

• Amplitudes of nonradial modes in classical Cepheids are very low, on aver- 
age almost an order of magnitude lower than in RR Lyr stars. 

• In RR Lyr variables nonradial modes are detected both in fundamental (RRab) 
and in first overtone pulsators (RRc). In fact, in all studied stellar sys- 
tems they are found in the RRab stars either equally frequently (Soszynski 
et al. 2002) or much more frequently (Soszynski et al. 2003; Mizerski 2003). 
In sharp contrast, in classical Cepheids nonradial modes are detected only 
in the first overtone pulsators. The difference between fundamental mode 
Cepheids and RR Lyr stars is even more striking, considering that for over- 
tone pulsators incidence rates of nonradial modes are roughly the same: 9.1% 
in FO Cepheids and 9.6 - 12.1 % in RRc stars (Nagy & Kovacs 2006, Mi- 
zerski 2003). 

• When two secondary frequencies are found in an RR Lyr star, they usually 
form together with the primary frequency an equally spaced triplet, centered 
on the primary (radial) peak. Such equidistant triplets are never observed in 
Cepheids. 

• In RR Lyr stars nonradial modes are detected only in close vicinity of the 
primary pulsation mode. This is also the case for most of the FO-v Cepheids, 
but in several variables a high frequency mode with Pv/P\ — 0.60 — 0.64 is 
found. Such high frequency modes were not detected in RR Lyr pulsators. 

5.6. Can secondary frequencies result from blending ? 

In Section 2 we describe how to identify spurious secondary frequencies, which 
have appeared through a contamination of DIA photometry with light of neigh- 
bouring variables. This procedure cannot identify cases of real blending, resulting 
from the angular resolution hmit of the observations. This raises a legitimite ques- 
tion, whether secondary periodicities discussed in this paper originate in Cepheids 
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or whether they are introduced by blending of Cepheids with other variable stars. 
Although the latter possibihty cannot be ruled out in any individual case, we are 
convinced that majority of nonradial modes detected in FO-v Cepheids cannot be 
explained by blending. Our arguments are statistical. First, observed distribution 
of frequency differences, A/ = /v — /i , is not random (Fig. 7). There is no reason 
why blending should introduce frequencies only in the vicinity of a radial mode or 
at any particular narrow range of period ratios. Second, we do not find nonradial 
modes in Cepheids, which pulsate in the fundamental mode. There is no reason 
why blending should affect stars pulsating in one mode, but not in the other mode. 
Third, blending should be more substantial (i.e. relative flux contribution should be 
larger) in case of fainter Cepheids (Mochejska 2002). Therefore, secondary perio- 
dicities introduced by blending should be found preferentially in short period pul- 
sators (where their amplitudes should be higher). The statistics of FO-V Cepheids 
shows just the opposite trend (Fig. 10). Finally, probability of blending should in- 
crease proportionally to the density of stars in the field (Kiss & Bedding 2005). 
This is not the case for the FO-v Cepheids. In 8 densest LMC fields (SC2-SC9; 
189 FO Cepheids in total) and in 8 least dense LMC fields (SC12, SC14, SC15, 
SC17-SC21; 157 FO Cepheids in total) the incidence rates of FO-V pulsations are 
10.6 lb 2.2% and 8.9 it 2.3%, respectively. The two values are statistically the 
same within 0.6CT, despite two-fold difference in average number of stars per field. 
From all these evidences we conclude that nonradial modes detected in the FO-V 
Cepheids do not result from blending, but are intrinsic to Cepheids themselves. 



6. Nonradial modes in FU/FO double-mode Cepheids 



OGLE-II catalog of double-mode LMC Cepheids (Soszyhski et al. 2000) lists 
19 FU/FO pulsators. Discovery of four new members (Section 3) brings the total 
number of Cepheids in this class to 23. We found nonradial modes in three of 
them. These are the first detections of nonradial modes in the FU/FO double-mode 
Cepheids. In the following, we will call these variables FU/FO-v Cepheids. The 
stars are listed in Table 6. In the fifth column of the Table we give the frequency 
difference between the nonradial mode and the first radial overtone. A/ = /v — 
/i . The prewhitened power spectra of the FU/FO-V variables are displayed in 
Fig. 1 1 . In the first two Cepheids, the secondary mode appears very close to the first 
overtone radial mode. The values of A/ are very similar to those observed in the 
FO-V Cepheids. In the third star, the secondary mode is found at high frequency, 
with Pv/^i = 0.6229. This is the same puzzling period ratio, which has been 
found in several FO-v pulsators. Clearly, nonradial modes excited in the FU/FO- 
V Cepheids are somehow related to the first radial overtone and their frequencies 
seem to be drawn from the same distribution as in the case of the FO-v Cepheids. 
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Table6 

FU/FO-V double-mode Cepheids in OGLE-H LMC sample 



OGLE ID 


Po [day] 


Pi [day] 


Pv [day] 


A/[c/d] 


Pv/Pl 


Av/Ai 


SC6-6 

SC15-69667 

SC17-126402 


4.841143 
1.904133 
3.209709 


3.453214 
1.376951 
2.294942 


3.581206 
1.517867 
1.429530 


-0.01035 
-0.06742 
0.26379 


1.03706 
1.10234 
0.62290 


0.057 
0.052 
0.040 



7. Blazhko effect in FO/SO double-mode Cepheids 

Including one new variable discovered in the current paper (Section 3), there are 
56 FO/SO double-mode Cepheids in the OGLE-11 catalog (Soszynski et al. 2000). 
In one-third of these stars we detected residual signal close to one or both primary 
(radial) pulsation frequencies. In all cases, however, the secondary frequencies 
were not resolved from the primary ones. Therefore, the frequency analysis of all 
FO/SO pulsators was repeated with MACHO b photometry (Allsman & Axelrod 
2001), which provided ~2.5 times better frequency resolution than OGLE-II. The 
MACHO data was cleaned in a manner described in Section 5.2 and then analyzed 
in a standard way. At this stage the OGLE-II sample was supplemented by 51 
additional LMC FO/SO double-mode Cepheids discovered by MACHO outside the 
area covered by the OGLE-II survey (Alcock et al. 1999, 2003). Our final FO/SO 
Cepheid sample consisted of 107 objects. 

With the MACHO data, residual power close to the primary pulsation frequen- 
cies was detected in 37 FO/SO double-mode Cepheids (35% of the sample). In 
20 of these stars, which constitute 19% of the sample, we were able to resolve 
this power into individual frequencies. Two frequencies are considered resolved 
if they differ by A/ > 0.000715 c/d {i.e. 1 /A/ < 1400 days). All resolved FO/SO 
Cepheids are listed in Table 7. 

7.1. Frequency domain 

Resolved FO/SO double-mode Cepheids display a very characteristic frequen- 
cy pattern. In most cases, in the vicinity of a radial frequency we detect two 
secondary peaks. They are located on the opposite sides of the primary (radial) 
peak and together with the primary peak they form an equally spaced frequency 
triplet. Such a structure is usually found around both radial modes. An example 
of this pattern is shown in Fig. 12. In several cases instead of triplets we find only 
doublets. This happens in objects, in which secondary frequencies are detected 
with the lowest signal-to-noise ratio. Therefore, it is most likely that these dou- 
blets are in fact parts of triplets, with the missing third component hidden in the 
noise. When signal-to-noise ratio is high we encounter an opposite situation. In 
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two FO/SO Cepheids we find equally spaced frequency quintuplets centered on the 
second overtone. Consecutive prewhitening reveals incomplete quintuplets in six 
other stars, including three variables where they apper around the first overtone. In 
Fig. 13 we display representative examples of triplets and quintuplets observed in 
the FO/SO pulsators. 

The secondary peaks are always very small, with amplitudes of 15mmag on 
average and exceeding 40 mmag only in one star. The side peaks of the multiplets 
are usually not equal. For the first overtone, the low frequency components of 
triplets/quintuplets are always higher than the high frequency ones. This is also 
true for 74% of multiplets centered on the second overtone. 

The separations of components in the two triplets/quintuplets, A/i and A/2, 
are also very small and never exceed 0.00143 c/d. This corresponds to the beat (or 
modulation) period of Pb = I /A/ > 700 day. The lower limit for A/ (upper limit 
for Pb) is currently given by resolution of the data. The physical upper limit for the 
beat (modulation) period Pg is unknown. 

When present around both radial modes, the two multiplets have nearly the 
same frequency spacings. Specifically, the difference between the two spacings, 
IA/2— A/i| is always below 8 x 10^^ c/d. According to Monte Carlo simulations 
of Alcock et al. (2000, 2003), for MACHO data such a difference is well within 
deviations expected from observational noise. In other words, triplets/quintuplets 
around both radial modes of a FO/SO Cepheid have the same frequency separation 
within accuracy of the data. 

Frequency separation A/ can take any value below 0.00143 c/d, but some va- 
lues are more likely than others. This is clearly visible in Fig. 14. Although statis- 
tics is rather small, the distribution of A/ shows a bimodal shape with pronounced 
maxima at ~ 0.00095 c/d and ~ 0.00125 c/d. These correspond to preferred beat 
(modulation) periods of Pb ~ 1050 days and Pb ^ 800 days. 



Table? 

Blazhko FO/SO double-mode Cepheids in combined OGLE-n + MACHO LMC sample 























^2 


'^2~ 




OGT F TD 


p, 
^ 1 


Pi 


Afr 
1 




Pb 


At 


At 


^1 




At 
^2 


At+ 
^2 






[day] 


[day] 


[c/d] 


[c/d] 


[day] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 




SC1^845 


0.951981 


0.766013 


0.00127 


0.00125 


794 


0.165 


0.021 




0.046 


0.021 




















0.011 






0.010 






SC 1-285275 


0.856624 


0.689234 


0.00112 


0.00109 


902 


0.177 


0.031 




0.040 


0.022 


0.008 


















0.014 






0.013 






SC 1-335559 


0.749802 


0.603647 


0.00121 


0.00128 


789 


0.227 


0.019 




0.062 


0.023 






SC2-55596 


0.932530 


0.751387 


0.00130 




768 


0.146 


0.007 




0.039 






a 
















0.007 


— 




— 


— 




SC6-142093 


0.896289 


0.722058 


0.00091 


0.00091 


1103 


0.154 


0.028 




0.043 


0.017 


0.012 


















0.013 






0.017 


0.006 




SC6-267410 


0.888530 


0.716773 




0.00117 


857 


0.120 






0.036 




























0.033 






SC8-10158 


0.689991 


0.555699 


0.00086 


0.00094 


1082 


0.175 


0.011 




0.031 


0.019 




c, si 






















0.012 






SC 11-233290 


1.217532 


0.978435 


0.00099 


0.00104 


991 


0.186 


0.019 




0.043 


0.017 


0.006 


















0.010 






0.010 






SC15-16385 


0.990417 


0.795736 


0.00084 


0.00088 


1185 


0.258 


0.017 




0.049 


0.020 




















0.016 












SC20-1 12788 


0.737740 


0.594287 


0.00072 


0.00072 


1390 


0.165 


0.062 


0.019 


0.015 


0.016 




















0.019 






0.013 






SC20-138333 


0.859788 


0.692171 


0.00127 


0.00126 


793 


0.189 


0.013 




0.060 


0.018 




















0.008 






0.010 







Table? 
Concluded 

















AT 






A2 


Ar 




MACHO ID 


Pi 


Pi 




A/2 


Pb 


Ai 


At 

1 


A++ 
1 




At 

2 


A++ 

2 


Remarks 




[day] 


[day] 


[c/d] 


[c/d] 


[day] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 




2.4909.67 


1.084120 


0.869989 


0.00099 


0.00098 


1012 


0.216 


0.012 

0.010 


— 


0.055 


0.013 

0.009 


— 


a,b 


13.5835.55 


0.898730 


0.722801 


0.00094 


0.00095 


1061 


0.244 


0.037 
0.021 


0.010 


0.040 


0.026 
0.016 


0.012 


a 


14.9585.48 


0.935802 


0.752812 


0.00091 


0.00092 


1098 


0.139 


0.024 
0.008 




0.036 


0.025 
0.010 


0.012 




17.2463.49 


0.762933 


0.614047 


0.00094 


0.00093 


1072 


0.235 


0.018 
0.014 




0.051 


0.013 
0.015 




d 


18.2239.43 


1.364186 


1.093324 


0.00142 


0.00142 


706 


0.230 


0.035 
0.014 


0.010 


0.063 


0.017 
0.015 




e 


22.5230.61 


0.633077 


0.510128 




0.00124 


805 


0.221 






0.045 


0.016 
0.014 






23.2934.45 


0.734350 


0.591775 


0.00129 


0.00125 


792 


0.207 


0.023 




0.055 


0.031 
0.026 


0.015 
0.011 


a 


23.3184.74 


0.841151 


0.677788 


0.00089 


0.00093 


1099 


0.160 


0.015 




0.048 


0.011 
0.013 






80.7080.2618 


0.715873 


0.577961 


0.00109 


0.00110 


914 


0.168 


0.010 




0.055 


0.006 
0.009 




d 



Remarks: a) Unresolved residual power at /i ; b) Unresolved residual power at /2; c) /j marginally detected; d) marginally 
detected; e) Additional peak at A/2 = -0.00076 [c/d]; si) Same as SC9-286128. 
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7.2. Time domain 

The observed triplet/quintuplet frequency pattern is a Fourier representation of 
a periodic modulation of both radial modes, with the common period Pg = 1/A/. 
In Fig. 15 we display this modulation for one of the FO/SO Cepheids. The plot 
has been constructed by dividing the data into 10 subsets, each covering 10% of 
the modulation phases, and then fitting radial frequencies and their linear combi- 
nations (Eq.(2)) to each subset separately. Fig. 15 shows that both amplitudes and 
phases of the radial modes undergo modulation. The amplitude variability is much 
stronger for the second than for the first overtone (±66% v*. ±17% change). Mini- 
mum amplitude of one mode coincides with maximum amplitude of the other mode. 
Maximum phase of the first overtone occurs just before maximum of its amplitude. 
For the second overtone, maximum phase coincides with minimum amplitude. We 
note, that the amplitude of the second overtone is modulated in a distinctively non- 
sinusoidal way. This explains presence of quintuplets in the frequency spectrum of 
this mode. 

Characteristic pattern of periodic modulation displayed in Fig. 15 is common 
to all FO/SO Cepheids of Table 7. The phenomenon is strikingly similar to the 
Blazhko modulation observed in the RR Lyrae stars (see e.g. Kurtz et al. 2000). 
Therefore, variables of Table 7 will be called Blazhko double-mode Cepheids, or 
FO/SO-BL Cepheids. 

7.3. FO/SO Cepheids with unresolved residual power 

In 17 FO/SO double-mode Cepheids (16% of the sample) we find after prewhi- 
tening a significant remnant power unresolved from the primary pulsation frequen- 
cies. This is a signature of slow changes of amplitudes and/or phase of the modes, 
not resolved within the length of the data. We call these stars FO/SO-PC Cepheids. 
Their list is given in Table 8. 

The amplitude and phase variability of each FO/SO-PC Cepheid was examined 
with method described in Section 5.2: we divided the data into 10 subsets, each 
spanning 10% of the total timebase, and then fitted Eq.(2) separately to each sub- 
set. In some stars we found only changes of pulsation phases, but in 10 variables 
we also detected clear changes of both amplitudes. In Fig. 16 we show typical ex- 
ample of such behaviour. The pattern of amplitude and phase variability is very 
much the same as that displayed in Fig. 15. In particular, the amphtudes of the two 
radial modes vary in opposite directions. This similarity suggests, that most of the 
FO/SO-PC Cepheids experience the same type of periodic modulation as Blazhko 
FO/SO Cepheids listed in Table 7, but on timescales longer than the current data. 

7.4. Incidence rate of Blazhko effect in FO/SO double-mode Cepheids 

Periodic amplitude and phase modulation in FO/SO double-mode Cepheids is 
by no mean a rare phenomenon. According to our analysis, its incidence rate is at 
least 18.7 ± 3.8%. The true rate can be perhaps even as high as 35%, if FO/SO-PC 
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Tables 

FO/SO-PC double-mode Cepheids in combined OGLE-H + MACHO LMC sample 



OGLE /MACHO ID 


Pi [day] 


P2 [day] Ai [mag] 


A2 [mag] 


Remarks 


SC4-176400 


1.108872 


0.895139 


0.089 


0.023 


ab 


SC4-220148 


0.740833 


0.595133 


0.222 


0.063 


ab 


SC7-120511 


1.251130 


1.003698 


0.133 


0.028 


ab 


SClO-204083 


0.526261 


0.422879 


0.245 


0.030 


b 


SCI 6-266808 


1.352920 


1.081016 


0.167 


0.025 


b 


SC17-186042 


0.610998 


0.491213 


0.179 


0.039 


ab 


SC20-1 88572 


1.015249 


0.814350 


0.167 


0.032 


ab 


SC21-12012 


1.341463 


1.074876 


0.181 


0.057 


ab 


6.6934.67 


0.920169 


0.740010 


0.181 


0.050 


ab 


11.9348.78 


0.737799 


0.594411 


0.200 


0.041 


a, si 


14.9098.35 


0.727783 


0.586544 


0.208 


0.032 


ab 


15.10428.60 


0.651848 


0.525346 


0.240 


0.044 


a 


23.3061.82 


0.595737 


0.480997 


0.166 


0.038 


ab 


24.2853.69 


0.708844 


0.571663 


0.269 


0.030 


b 


24.2855.80 


0.593778 


0.478320 


0.178 


0.020 


b 


47.2127.102 


0.578454 


0.466643 


0.227 


0.024 


ab 


80.7079.62 


1.347888 


1.076366 


0.177 


0.042 


b 



Remarks: a) Residual unresolved power at /i; b) Residual unresolved power 
at /i; si) Same as 14.9348.3191. 



variables are confirmed to undergo Blazhko modulation with very long periods, as 
we suspect. 

In Fig. 17 we present the Petersen diagram for our sample of FO/SO double- 
mode Cepheids of the LMC. Blazhlco FO/SO pulsators and FO/SO-PC pulsators 
are marlced with filled and open circles, respectively. Blazhko FO/SO variables 
are not detected at all for Pi < 0.6 day, but for longer periods they intermingle with 
"normal" double-mode Cepheids. The same can be said about FO/SO-PC variables, 
except that the short period cutoff is at Pi = 0.5 day. Apart from amplitude and 
phase variability, Blazhko pulsators and PC pulsators do not differ from the rest of 
FO/SO double-mode Cepheids. 

Fig. 17 shows, however, that the incidence rate of Blazhko effect might depend 
on pulsation period. To address this point we display in Fig. 18 the histogram of 
first overtone periods for our sample of FO/SO Cepheids. The incidence rate of 
Blazhko effect peaks at 0.8 day < Pi < 1.0 day, where it reaches 48%. It sharply 
dechnes below 15% for both longer and shorted periods. While the decline of the 
incidence rate towards shorted periods (i.e. lower luminosities) can be attributed to 
the selection effects, the decline towards longer periods must be real. 
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7.5. What Causes the Modulation ? 

Any model of Blazhko FO/SO double-mode Cepheids has to explain two most 
basic properties of these stars: 

• both radial modes are modulated with the same period. 

• the amphtude variations of the two modes are anticorrelated: maximum of 
one amphtude coincides with minimum of the other. 

Many different ideas have been put forward to explain Blazhko modulation in 
monoperiodic RR Lyrae stars (see Stothers 2006 for most recent summary). Two 
primary contenders, which are most popular nowadays, are the oblique magnetic 
pulsator model (Shibahashi 1995, 2000) and the 1: 1 resonance model (Nowakowski 
& Dziembowski 2001). Another, very different scenario has been proposed re- 
cently by Stothers (2006). We will argue, that all these models fail to account for 
modulation observed in the FO/SO double-mode Cepheids. 

Shibahashi's obhque magnetic pulsator model assumes presence of a dipole 
magnetic field inclined to the rotation axis of the star. The field introduces quadru- 
pole distortion to the radial mode. As the star rotates, the distortion is viewed from 
different aspect angles, which leads to variation of an apparent pulsation ampli- 
tude. Somewhat similar idea has been considered by Kovacs (1995; see also Balazs 
1959). In his scenario the distortion is caused by an £ = 1 nonradial mode, which 
is aligned with a magnetic axis of the star. The nonradial mode is excited through a 
1 : 1 resonance with the radial pulsation. Consequently, frequencies of both modes 
are exactly synchronized. Both in Shibahashi's and in Kovacs's models physical 
amplitudes of pulsation are constant. Modulation of the observed amplitudes is 
caused by geometrical effect of rotation. The models naturally explain why both 
radial modes in FO/SO Cepheid are modulated with the same period. However, in 
the above scenario both modes should reach maximum amphtudes simultaneously, 
because they are both distorted in the same way. This is not what we observe. Thus, 
the Blazhko effect in FO/SO double-mode Cepheids cannot be reproduced by the 
oblique magnetic pulsator. 

The model proposed by Nowakowski & Dziembowski (2001, see also Kovacs 
1995) assumes a 1:1 resonant couphng of the radial mode to a pair of nonradial 
modes of C = \ ,m = ±.\ . Such a mechanism generates a triplet of equally spaced 
frequencies. The physical amplitudes of the modes are constant, but beating of 
their equidistant frequencies leads to a slow periodic modulation of apparent am- 
phtude and phase of pulsation. The 1 : 1 resonance model naturally explains triplets 
observed in frequency spectra of Blazhko FO/SO Cepheids. Quintuplets can also 
be easily understood, by generahzing the model to couphng with modes of £ = 2. 
However, in the scenario of Nowakowski & Dziembowski (2001) modulation of 
each radial mode is an independent process. Therefore, there is no reason why both 
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radial modes should be modulated with exactly the same period (although modu- 
lation periods shouldn't be very different). There is also no reason why amplitude 
variations of the two modes should be in any particular phase relation to each other, 
as is observed. Finally, the 1 : 1 resonance model does not explain why in vast ma- 
jority of cases modulation is observed either for both radial modes or for none of 
them. 

The new model (or rather an idea) of Stothers (2006) fails to explain the ob- 
servations of Blazhko FO/SO Cepheids, as well. In his scenario it is assumed that 
turbulent convection in the stellar envelope is cyclically weakened and strengthened 
by a transient magnetic field generated by dynamo mechanism. As convection has 
a strong effect on the amplitudes of pulsation (e.g. Feuchtinger 1999), the latter 
become modulated, too. In this picture we would expect both modes in a FO/SO 
Cepheid to reach maximum amplitudes at the same time, because both amplitudes 
vary in response to a common factor - a changing strength of convection. This 
prediction is in obvious conflict with observations. 

At this stage, the mechanism causing amplitude and phase modulation of FO/SO 
double-mode Cepheids remains unknown. Nevertheless, the available observations 
already provide some important clues. The fact that the two modes are always 
modulated with the same period proves that their behaviour is not independent. 
Both modes must be part of the same dynamical system. The fact that high ampli- 
tude of one mode always coincides with low amplitude of the other mode strongly 
suggests that energy transfer between the two modes is involved. Thus, available 
evidences point towards some form of mode interaction, which causes a periodic 
amplitude and phase modulation of the modes. Such a modulation (limit cycle) 
can be induced by a resonant coupling of a radial mode with another radial or 
nonradial mode (Moskalik 1986; Van Hoolst 1995). Because of the nonlinear fre- 
quency synchronization, the resonant mode will not appear as an independent peak 
in the power spectrum. It is enough that only one of the observed radial modes is 
modulated by the resonant coupling. This modulation will be carried over to the 
other radial mode through the so-called cross-saturation effect. Speaking in physi- 
cal term, the two radial modes compete for the same driving (k mechanism in the 
He"*" ionization zone). When the amplitude of one mode is suppressed by whatever 
mechanism (e.g by a resonance), it allows the other mode to grow. This explains 
in a very simple and natural way why the two radial modes are modulated with the 
common period and why their amplitudes are always anticorrelated. 

8. Conclusions 

Taking advantage of a large and homogenous photometric dataset collected 
during OGLE-II project, we have performed a systematic frequency analysis of 
classical Cepheids identified in the Large Magellanic Cloud (Udalski et al. 1999b; 
Soszyiiski et al. 2000). This study has been aimed at finding multiperiodic and 
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nonstationary variables among Pop. I Cepheids and at establishing reliable statis- 
tics of different forms of their pulsational behaviour. The main body of our survey 
has been conducted with OGLE-II photometry, but in cases when frequency reso- 
lution turned out to be insufficient, we have resorted to photometry collected by the 
MACHO experiment. 

We have discovered two triple-mode Cepheids, which pulsate with three lowest 
radial overtones simultaneously excited. Triple-mode radial pulsators have been 
known before, but only among High Amplitude 8 Scuti stars (Kovacs & Buchler 
1994; Antipin 1997b; Handler, Pikall & Diethelm 1998; Wils et al. 2008). This is 
the first detection of such pulsators among Cepheids. Three additional triple-mode 
Cepheids have been identified in the LMC by Soszynski et al. (2008a). These rare 
variables are a very valuable trophy, since their seismological analysis can strongly 
constrain the stellar evolution theory (Moskalik & Dziembowski 2005). 

In 19% of FO/SO double-mode Cepheids we have detected periodic variability 
of amplitudes and phases of the two radial modes. Another 16% of FO/SO pul- 
sators is suspected of undergoing the same type of modulation, but on timescales 
longer than our data. This phenomenon is analogous to the Blazhko modulation 
in the RR Lyr stars, which was first discovered a century ago (Blazhko 1907). In 
Blazhko FO/SO Cepheids both modes are modulated with a common period, al- 
ways in excess of 700 days. The amplitudes of the modes vary in opposite phases, 
that is a maximum amphtude of one mode coincides with minimum amplitude of 
the other. In frequency domain, the modulation manifests itself as splitting of each 
radial mode into a triplet or sometimes a quintuplet of equidistant peaks, spaced by 
frequency of the modulation. 

The discovery of modulated FO/SO Cepheids shows that Blazhko effect and 
double-mode pulsations are not mutually exclusive. These stars will play a very 
important role in selecting a proper model to explain the Blazhko phenomenon. 
Observations of two modulated modes in one star will put any proposed scenario 
to a very stringent test. Two most popular models, the obUque magnetic pulsator 
model (Shibahashi 2000) and the 1:1 resonance model (Nowakowski & Dziem- 
bowski 2001) have already failed this test and should be ruled out, at least in the 
case of FO/SO Cepheids. 

Perhaps the most exciting result of our survey is detection of nonradial modes 
in classical Cepheids. Such modes have been found in 42 overtone pulsators (which 
constitutes 9% of the sample) and in three FU/FO double-mode pulsators. They 
have not been found in Cepheids pulsating in the fundamental mode. 

We find two different types of nonradial modes in Cepheids. In most cases these 
modes are detected very close to the first radial overtone, but in several stars they 
appear at considerably higher frequencies. In the later case, the measured period 
ratio falls in a narrow range of 0.60 — 0.64, which places the nonradial mode in the 
close vicinity of the (unobserved) fourth radial overtone. The observed frequency 
pattern cannot be explained by any form of amphtude or phase modulation, it is 
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also incompatible with excitation of radial modes. This makes interpretation in 
terms of nonradial modes the only possibility. 

Detection of nonradial modes in classical Cepheids has been claimed before 
by Kovtyukh et al. (2003), based on observations of bumps in spectral line profiles 
of four galactic Cepheids. Such structures can be indicative of nonradial oscilla- 
tions {e.g. Telting 2003). However, the authors have not shown that the observed 
line profiles vary in a periodic way, or that they vary with their own period differ- 
ent from that of the radial mode. This is an important test, because bumps in the 
spectral lines can be caused by mechanisms other than nonradial pulsations {e.g. by 
shocks). Therefore, observations presented by Kovtyukh et al. (2003) are in our 
opinion inconclusive. The results of frequency analysis discussed in the current 
paper provide the first solid evidence of nonradial modes excitation in classical 
Cepheids. 

Discovery of nonradial modes in Cepheids poses a major challenge to the stel- 
lar pulsation theory. In order to be photometrically detectable, nonradial modes 
have to be of low spherical degree, £ (Dziembowski 1977). However, linear pul- 
sation calculations predict, that in classical Cepheids all modes of I < A should 
be heavily damped (Osaki 1977). This conclusion has recently been confirmed by 
new calculations of Mulet-Marquis et al. (2007), who have found that all nonradial 
modes of £ < 5 should be damped. We note, that this is a very different situation 
from that encountered in the RR Lyr models, where many modes of ^ = 1 , 2 are 
linearly excited (Dziembowski & Cassisi 1999). 

When writing of this manuscript was almost completed, we learnt about a new 
preprint on LMC Cepheids by Soszyriski et al. (2008b). The authors have analyzed 
data collected during the third phase of the OGLE project (OGLE-HI). Among 
many interesting results, they have also detected secondary periodicities in all sub- 
classes of Pop. I Cepheids. This provides an important confirmation of our results 
with independent data. The incidence rates given by Soszynski et al. are very simi- 
lar to ours in case of the double-mode pulsators, but significantly higher than ours 
in case of the FO Cepheids. This is most likely due to much longer time span 
of their photometry and higher number of measurements. Interestingly, Soszynski 
et al. has detected secondary periodicities also in the FU Cepheids, which we have 
not. A detailed comparison of incidence rates derived by us and by Soszynski 
et al. is beyond the scope of the present paper. We only note, that such compari- 
son requires some caution, because the latter authors do not discriminate between 
secondary periodicities resolved from the primary frequency and those, which are 
not resolved. 
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Fig. 1. Petersen diagram for FU/FO double-mode Cepheids in the Galaxy (open squares), LMC 
(asterisks) and SMC (open circles). New LMC double-mode Cepheids are displayed with filled 
circles. Dotted lines represent the best linear fits to the period ratios observed in each population. 
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Fig. 2. Petersen diagram for FO/SO double-mode Cepheids in the LMC. New double-mode Cepheid 
SC20-112788 displayed with a filled circle. 
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Fig. 3. Prewhitening sequence for triple-mode Cepheid SC5-338399. Upper panel displays power 
spectrum of original data. Lower panels show power spectra after removing consecutive frequencies. 
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Fig. 4. Wi index P — L relations for LMC Cepheids. Upper and lower sequences correspond to 
first overtone and fundamental mode Cepheids, respectively. First overtone periods of triple-mode 
Cepheids SC3~360128 and SC5-338399 are displayed with black symbols. 
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Fig. 5. Prewhitened power spectra of sample of FO-V Cepheids with small |A/| . Frequencies of 
removed radial modes indicated by dashed Unes. 
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Fig. 6. Prewhitened power spectra of sample of FO-v Cepheids with Pv/Pi = 0.60 — 0.64. Frequen- 
cies of removed radial modes indicated by dashed lines. 
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Fig. 7. Distribution of frequency differences Af = fy—fi for LMC FO-V Cepheids and for LMC 
RRc-BL stars (from Nagy & Kovdcs 2006). For RRc-BL variables with frequency triplets, only 
higher side peaks are taken into account. Shaded area marks variables with Py/Pi =0.60 — 0.64. 
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Fig. 9. W; index P — L diagram for LMC Cepiieids. Dominant (radial) periods of FO-V Cepheids 
are displayed with black symbols. 
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Fig. 10. Incidence rate of LMC FO-v Cepheids vs. first overtone period. 
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Fig. 11. Power spectra of FU/FO-v double-mode Cepheids after prewhitening with frequencies of 
the two radial modes and their linear combinations. Removed radial frequencies indicated by dashed 
lines. 
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Fig. 12. Prewhitened power spectrum of Blazhko FO/SO double-mode Cepheid SCl-44845. Fre- 
quencies of removed radial modes indicated by dashed Unes. Lower panels display the fine structure 
around the radial modes. 
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Fig. 13. Prewhitened power spectra of sample of Blazhko FO/SO double-mode Cepheids. Frequen- 
cies of removed radial modes indicated by dashed Unes. 
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Fig. 14. Distribution of modulation frequencies, Af, for Blazhko FO/SO double-mode Cepheids. 
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Fig. 15. Periodic modulation of amplitudes and phases of radial modes in FO/SO double-mode 
Cepheid SCl-285275. First and second overtones are displayed with filled and open circles, respec- 
tively. The error bars are ±lo. Modulation period is = 891.6 day. 
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Fig. 16. Amplitude and phase variations in FO/SO-PC double-mode Cepheid SC17-186042. Sym- 
bols are the same as in Fig. 15. 
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Fig. 17. Petersen diagram for FO/SO double-mode Cepheids in the LMC. Cepheids with Blazhko 
effect and FO/SO-PC variables displayed with filled and open circles, respectively. Remaining pul- 
sators plotted with asterisks. 




Fig. 18. Distribution of first overtone periods of FO/SO double-mode Cepheids in the LMC. Black 
histogram represents Blazhko pulsators. Shaded area marks FO/SO-PC variables. 



